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ABSTRACT 

The transient neutron star low-mass X-ray binary EXO 1745-248, located in the 
globular cluster Terzan 5, was detected during its quiescent state with Chandra 
in 2003. The source displayed a 0.5-10 keV quiescent X-ray luminosity of L q ~ 
10 33 (D/5.5 kpc) 2 erg s _1 , which was completely dominated by hard non-thermal 
emission. This is at odds with other non-pulsating neutron stars that typically show 
detectable soft thermal emission at such quiescent luminosities. Here we use three 
additional Chandra observations, performed in 2009 and 2011, to further study the 
quiescent properties of EXO 1745-248. We find that the powerlaw intensity varies 
considerably up to a factor of ~ 3 within hours and by about one order of magnitude 
between the different epochs. We discuss the implications of the observed change in 
quiescent flux for the interpretation of the hard powerlaw emission. We constrain the 
neutron star surface temperature as seen by a distant observer to kT°° < 42 eV and 
the thermal bolometric luminosity to L qj boi ^ 7 x 10 31 (D/b.5 kpc) 2 erg s^ 1 . This 
confirms that EXO 1745-248 harbours a relatively cold neutron star and suggests 
that, for example, enhanced cooling mechanisms are operating in the stellar core or 
that the binary on average resides in quiescence for hundreds of years. 
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- X-rays: individual (EXO 1745-248) - accretion, accretion discs 



1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are binary star systems 
that are composed of a compact primary (a neutron star 
or a black hole) that accretes the outer gaseous layers of 
a (sub-) solar-mass companion star. Neutron star primaries 
can reveal themselves by showing coherent X-ray pulsations 
when the accretion flow is funnelled onto the magnetic poles 
of the neutron star, or by displaying type-I X-ray bursts 
caused by unstable thermonuclear burning of the accreted 
matter on the surface of the neutron star. 

Many neutron star LMXBs are transient and alternate 
active accretion outbursts with long episodes of quiescence. 
Whereas outbursts typically last a only few weeks and gen- 
erate O^UTkeV X^ luminosities of L x ~ 10 36 " 38 erg s" 1 
(e.g., IChen et al.lll997h . the quiescent phase may extend to 
many decades and is characterised by a much lower qui- 
escent luminosity of L q ~ 10 31 ~ 33 erg s _1 (0.5-10 keV; 
iHeinke et al.ll2009h . The cause of th e quiescent X-ra y emis- 
sion has been subject of debate (e.g.. lCampanall2003h . Spec- 
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trally, one can distinguish a soft thermal component at en- 
ergies below ~ 2 keV, and a hard non-thermal emission tail 
that dominates the quiescent X-ray spectrum above ~ 2 — 3 
keV. 

The soft quiescent spectral component is generally in- 
terpreted as heat being radiated from the surface of the cool- 
ing neutron star ijBrown et al.lll99Sl : IColpi et all 1200 ll ). As 
such, the quiescent thermal emission provides a measure of 
the neutron star's interior temperature. However, low-level 
accretion onto the surface of the neutron star might a lso 
produce a thermal X-ray spectrum (|Zampieri et al.l ll995). 

The hard non-thermal emission is usually modelled as 
a simple power law with a photon index of V ~ 1 — 2, and 
its fractional contribution to th e total 0.5-10 keV fl ux may 
vary anywhere from — 100% (|jonker et al.ll2004l ). It has 
been attributed e.g., to shock emission generated in the in- 
teraction of a residual accretion flow with the magnetic field 
of the neutron star, or the re-activa tion of a radio pulsar 
in quiescence l|Campana et al.lll998l ). It remains an unan- 
swered question if, and how, the two spectral components 
are related. 
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Table 1. Historic activity from Terzan 5. 



Year Satellite(s) 


Lx 


-^X,peak 


i b Comments 


Ref. 




(erg s" 1 ) 


(erg s" 1 ) 


(days) 





1980 


Hakucho 


< 3 x 10 36 




> 16 


Type-I X-ray bursts 


1,2 


1984 


EXOSAT 


- 3 x 10 37 








3 


1990 


ROSAT (PSPC) 


~ 3 x 10 35 








1 


1991 


ROSAT (HRI) 


~ 3 x 10 34 








5 


2000 


RXTE, Chandra 


- 2 x 10 37 


~ 6 x 10 37 


~ 100 


Type-I X-ray bursts, dipping 


6,7,8 


2002 


RXTE 


~ 2 x 10 37 


~ 4 x 10 37 


~ 30 




8,9 


2010 


Integral, RXTE, Swift, Chandra 


~ 2 x 10 37 


~ 7 x 10 37 


~ 80 


Type-I X-ray bursts, 11-Hz pulsations 


10-15 


2011 


MAXI, RXTE, Swift, Chandra 


~ 3 x 10 36 


~ 9 x 10 36 


~ 20 


Superburst 


16-19 



Note. - Lx represents the (average) 0.5-10 keV luminosity and t j, reflects any available constraints on the outburst du- 
ration. We calculated Lx from fluxes reported in literature or using publicly available RXTE / ASM or MAXI data. For 
the conversion to the 0.5—10 keV energy band we use PIMMS (v. 4.4) by assuming a fiducial spectrum with Ng = 1 X 
10 22 cm -2 and T = 1.5, and w e adopt a distance of D = 5.5 kpc. When poss ible we also gi v e the 0.5 -10 keV peak 
lumin osity Lx.peak- Refer e nces: l= lMakishima et al.l lll98ll ). 2=llnoue et al.l lll984h 3=IWa.rwick et alj lll988l ). 4= IVerbunt et al.l 
l|l995h . 5=ljohnston et all lll995l) . 6 =jMarkwardt fc SwanM ll2000h. 7=lHeinke et al. |2003), 8=Altam ir ano e t al. in prepara- 
tion, 9= IWiinands et alj 2002h . 10 =lBordas et alj j201oF~ ll=|Papitto et al.l (l20lTh. 12 = jMiller et al l ll201ll). 13=lMotta et al.l 



J201lf). 14=lLmares etldl feoill). 15=lDegenaar fc Wiinandsl ll2011ah . 16= lAltamirano et alj ll2011ah . 17= lAltamirano et alj fcOllbl) , 



18= lMihara et al.l ||20hT k 19= |Poolev et alj j201lF 



1.1 EXO 1745-248 in Terzan 5 

Terzan 5 is a globular cluster located in the bulge of 
our Galaxy at an e stimated distance of D = 5.5 kpc 
l|Ortolani et a.1.1 [20071 ') PI This cluster has a particularly high 
stellar density, which is thought to cause a high rate of dy- 
namical star encou nters and the form ation of compact bi- 
nary systems (e.g.. IPoolev et al"1l2003l ) . Apart from a large 
number of millisecond radio pulsars (I Ransom et al.1 120051 ) 
and faint X-ray point sou rces that may be qu iescent LMXBs 
or cataclysmic variables (IHeinke et al.l 120061 ). Terzan 5 also 
contains two confirmed transiently accreting neutron stars. 

Already in 1980, a number of type-I X-ray bursts was 
detected from the direction of the cluster, which testified the 
prese nce of an active neutron star LMXB ( Mak ishima et all 
Il98ll ). Subsequent X -ray activity was observed from Terzan 
5 in 1984, 1990, 1991, 2000, 2002, 2010 and most recently in 
2011 (TablefTJand references therein). While the detections 
in 1984, 1990 and 1991 concern single observations and hence 
only snapshots of the activity, the most recent four outbursts 
have been covered by several X-ray instruments, e.g., by the 
ASM onboard RXTE (Fig. [[J. 

Multiple type-I X-ray bursts were detected during 
the 2000 and 2010 outbursts, while a superburst (i.e., a 
very energetic thermonuclear burst that lasts for several 
hours) was observed in 2011. This thus evidences active 
accretion onto a neutron star in 2000, 2010 and 2011 
iGallowav et al.ll200Sl;lLinares et alhoiltlMotta et al.ll201ll ; 
lAltamirano et al.ll2011bl ; iMihara et alj|201ll ). High~patial 
resolution Chandra observations have revealed that these 
three outbu rsts involved two di f ferent transient neutron s tar 
LMXBs (cf. IHeinke et al.ll2003l ; IPoolev et al.ll2010l . l201ll ). 

Historically, the source that was active in 2000 and 2011 
is referred to as EXO 1745-248, whereas the recently discov- 



ered 2010 transient has been dubbed IGR J1748O-24460 In 
addition to exhibiting type-I X-ra y bursts, the latter also 
displays H-Hz X- ray pulsations (jStrohmaver et all l2010l : 
IPapitto et al.|[201ll ) . The two transient neutron star LMXBs 
have a spatial separation of merely ~ 5" (see Fig. [2} ■ Hence, 
apart from Chandra, none of the present and past X-ray mis- 
sions provides sufficient angular resolution to spatially re- 
solve them. It is thus unclear which of the two transient neu- 
tron star LMXBs was causing the outbursts that occurred in 
1980, 1984, 1990, 1991 and 2002. Moreover, 50 distinct low- 
luminosity X-ray point sources have b een identified withi n 
the half-mass radius of Terzan 5 (0.83'; IHeinke et alj|2006r i. 
so the possibility of additional active transient sources caus- 
ing these historic outbursts cannot be ruled out. 

Terzan 5 has been observed with Chandra on several 
occasions when no bright transients were active (Table [2}. 
These observations can be addressed to investigate the two 
neutron star LMXBs during quiescence. The quiescent prop- 
erties of the 2010 transient, IGR J17480-2446, are the 
subjec t of a recent ongoing study |Degenaar fc Wiinandsl 
l2011al fbl lDegenaar et alj|201ll ). The quiescent emission spec- 
trum of the 2000/2011 transi ent, EXO 1745-248, was stud- 
ied by IWiinands etafl l|2005l ) using a Chandra observation 
obtained in 2003. At that time the source displayed a 0.5-10 
keV luminosity of L q ~ 10 33 erg s^ 1 , which i s not unusual 
for n eutron star LMXBs in quiescence (e.g., IHeinke et al.l 
2009). Unlike other non-pulsating neutron star LMXBs, 
however, the quiescent X-ray spectrum of EXO 1745-248 
was found to be completely dominated by a hard emis- 
sion component, with no indications of thermal radiation 
I Wiinands et alj|2005l ). Although a non-thermal component 
is often seen for n eutron star LMXBs in quiescence (e.g., 
Ijonker et al.|[20o3 ). typically thermal emission can be de- 
tected as well. 

The lack of detectable quiescent thermal emission 
places strong constraints on the properties of the neutron 



1 Throughout this work we assume a distance of D = 5.5 kpc for 2 In the study of lHeinke et al ] l|2006l l. EXO 1745-248 is denoted 

Terzan 5 and EXO 1745-248. as CX3 and IGR J17480-2446 as CX25. 
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Figure 1. /JXTS/ASM 5-d averaged lightcurve of Terzan 5, illustrating three different outbursts that were recorded in the past 15 years 
(1.5—12 keV). Arrows indicate Chandra observations that were performed when no bright X-ray transients were active (Tabled- We 
note that the most recent outburst that was detected in 2011 October— November (EXO 1745-248) was not covered by the RXTE / ASM. 



star core Jjonker et all I2004]. 120071: iTomsick et alJ 120051 : 
IWiinands et alj 120051 ; iHeinke et all 120091 '). These neutron 
stars must be cold, suggesting that their interior may be effi- 
ciently cooled through neutrino emissions (see also Sec. 14. l[l . 
Such neutron stars are expected to be relatively massive, 
because larger central densities are thought to lead to a 
higher rate of neutrino co oling (|Lattimer fc Prakashl |2003 ; 
lYakovlev fc Pethicld liool ). In this work, we discuss three 
additional Chandra observations of Terzan 5, carried out in 
2009 and 2011, to further investigate the quiescent proper- 
ties of EXO 1745-248. 



2 CHANDRA OBSERVATIONS OF TERZAN 5 

We use four Chandra observations of Terzan 5, spread be- 
tween 2003 and 2011, to study any possible spectral and 
temporal variability in the quiescent emission of EXO 1745- 
248. Details of the individual observations can be found in 
Table [2] and references therein. All four were carried out in 
the faint data mode, with the nominal frame time of 3.2 s 
and the target positioned on the S3 chip. Data reduction was 
carried out using the CIAO tools (v. 4.3) and following stan- 
dard Chandra analysis threadsjj The 2003 and 2009 data 
were reprocessed using the task ACIS_PROCESS_EVENTS 
to benefit from the most recent calibration. Episodes of high 
background were removed from the 2003 data, which re- 
sulted in a net exposure time of 31.2 ks. No background 
flares were present during the 2009 and 2011 observations, 
so the full exposure time was used in the analysis. 

EXO 1745-248 is clearly detected in the dense core of 
the cluster in 2003 and 2009 (Fig. [3 left). We used a circular 
region with a 1.5" radius to extract source events, and one 



with a radius of 40" positioned on a source-free part of the 
CCD ~ 1.4' west of the cluster core as a background refer- 
ence. We extracted count rates and lightcurves using the tool 
DMEXTRACT, while the meta-task SPECEXTRACT was 
used to obtain spectra and to generate the ancillary response 
files (arf) and redistribution matrix files (rmf). The spectral 
data was grouped into bins with a minimum of 15 photons 
using the tool GRPPHA and fitted using the software pack- 
age XSpec (v. 12.7). To take into account the interstellar 
neutral hydrogen absorption along the line of sight, we in- 
clude the PHABS model in all our spectral fits using the 
default XSpec abundances and cross-sections. Throughout 
this work we assume a distance of D = 5.5 kpc towards EXO 
1745-248 when converting X-ray fluxes into luminosities. All 
quoted errors correspond to 90% confidence levels. 

Whereas EXO 1745-248 was one of the brightest sources 
in Terzan 5 during the 2003 and 2009 observations, the 
source had considerably faded in 2011 (Fig. [5J). Although 
the source is not clearly visible in the 2011 data sets, 
there appears to be an excess of photons present at the 
source position. To verify this, we employed the wave let 
detection algorithm PWDETECT (|Damiani et al.lll997al lbl). 
which has been found to be particularly effective in de- 
tecting faint sources located close to brighter objects and 
hence it is a very useful tool for globular clusters (e.g., 
IHeinke et ail 12009 ). We performed standard PWDETECT 
runs on the ACIS-S3 chip with wavelet sizes varying from 
0.5" - 2.0" Q This exercise suggests that EXO 1745-248 
was weakly detected during both 2011 observations (see 
Sect. 13. 2p . We extracted count rates and spectra using the 
tools DMEXTRACT and SPECEXTRACT, respectively. 
We combined the 2011 spectral data using the ClAO-task 
COMBINE_SPECTRA to improve the statistics. 



3 http: / /cxc. harvard.edu/ciao/threads/index. html. 
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Figure 2. Chandra/ ACIS images obtained in 2009 and 2011 (0.5—5 keV). For the 2011 image, we summed the two individual exposures 
of February and April (Table [2}. The locations of the two known transient neutron star LMXBs are indicated by circles of 1.5" radii. 



Table 2. Chandra observations of EXO 1745-248 in quiescence. 



Date 


Obs ID 


£cxp 


Net count rate 


Ref. 


(yyyy-mm-dd) 




(ks) 


(counts s" 1 ) 




2003-07-13/14 


3798 


39.5 


7.1 X 10" 3 


1,2,3 


2009-07-14/15 


10059 


36.4 


1.1 X 10" 2 


3 


2011-02-17 


13225 


29.7 


4.0 x 10" 4 


4 


2011-04-29/30 


13252 


39.5 


3.8 x 10" 4 


5 



Note. - After correcting for background flares, the exposure 



time of the 2003 observation is t c 



31.2 ks. The quoted net 



source count rates are for the 0.5-10 ke V band. Observatio n 
de tails can be found in th e re ferences: l=[H einke et al.l ll2006h . 



Wijnands et al.l J2005l) . 3=lDegenaar fc Wiinanddj2011bF) , 
Degenaar*^T^Viinandsl l|2011aT l. S ^DegenaarlrtTliU l|201lTl . 



3 EXO 1745-248 IN QUIESCENCE 

3.1 The 2003/2009 data: spectra and variability 

The 2003 Chandra o bservation of Terzan 5 was used by 
IWiinands et~afl (|2005l ) to study the quiescent properties of 
EXO 1745-248. At that time the source was detected at 
a mean count rate of (0.71 ± 0.05) x 10~ 2 counts s" 1 and 
displaying significant variability in the light curve (0.5-10 
keV). The quiescent spectrum was found to be hard and ab- 
sent of a thermal emission component. The data could be 
best described by a single absorbed powerlaw model with 
a hydrogen column density Nh = ( 1.4 ± 0.5) x 10 22 cm" 2 
and a spectral index of V = 1.8=1= 0.5. IWiinands et ail (|2005l ) 
concluded that any possible thermal emission component 
contributed at most ~ 10% to the unabsorbed 0.5-10 keV 
flux of F x = 2.2tol x lcr13 er S cm ~ 2 s _1 . An upper limit 
on the neutron star temperature as seen by an observer at 
infinity of kT°° < 80 eV was inferred. We reprocessed and 
refitted the 2003 data in this work and obtained similar re- 
sults as found previously (see Table [3]). 

During the 2009 observation EXO 1745-248 was de- 
tected at an average count rate of (1.09 ± 0.06) x 
10" 2 counts s" 1 (0.5-10 keV). The count rate lightcurve, 
shown in Fig. [3l reveals considerable variation by a factor of 



~ 3 on a time scale of hours. We fitted the 2009 spectral data 
using a variety of models. Similar to the results for the 2003 
observation, we find that a thermal emission model does not 
provide a good description of the data. We first tried a neu- 
tron star atmosphere model, which is often used to describe 
the X-ray spectra of quies cent neutron star LMXBs. Us- 
ing the model NSATMOS (|Heinke et al.ll2006l ) we explored 
fits by either fixing several parameters or leaving these free 
to float. All trials resulted in a reduced chi-squared value 
xt > 6, hence no acceptable fit could be obtained. 

A simple blackbody can describe the data, but returns 
unrealistic model parameters: Nh < 0.3 x 10 22 cm -2 , a 
temperature kT = 1.4 ± 0.2 keV and an emitting radius 
R = 0.04±™ 4 km (for D = 5.5 kpc), yielding X l = 1-1 
for 21 degrees of freedom (dof). The obtained hydrogen col- 
umn density is well below that inferred for the 2000 out- 
burst and 2003 quiescent data of EXO 1745-248 (N H ~ 
1.5 x 10 22 cm- 2 : lHeinke et al.ll2003l : IWiinands et al.ll2005l ). 
as well as the average value found fo r the X-ray source s 
in Terzan 5 (N H ~ 1.9 x 10 22 cm" 2 : iHeinke et al.l [20061 ). 
Furthermore, the temperature and emitting radius would 
be highly unusual for a (non-pul sating) thermally-em itting 
neutron star in quiescence (e.g., iRutledge et al.l 1999). We 
thus conclude that this model does not correctly describe 
the data. 

Fitting the 2009 spectral data with a simple absorbed 
powerlaw yields N H = (1.1 ± 0.3) x 10 22 cm" 2 and T = 
1.3 ± 0.3 (x 2 = 1.3 for 22 dof). The resulting 0.5-10 keV 
unabsorbed flux is F x = (2.7 ± 0.2) x 10" 13 erg cm" 2 s" 1 
(see Table[3]). To constrain any thermal emission component, 
we add an NSATMOS model with a canonical neutron star 
mass and radius of Mns = 1-4 Mq and Rns = 10 km, a 
source distance of D = 5.5 kpc and a normalisation of unity. 
The only free fit parameter for this model component is the 
neutron star temperature. Addition of this model compo- 
nent does not improve the fit (xt = 1.3 for 21 dof). Since 
the spectral data do not require the inclusion of a thermal 
component, we consider the obtained neutron star temper- 
ature and thermal flux as upper limits to their true values. 
As such we find that the thermal component contributes 
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< 42% to the total 0.5-10 keV unabsorbed flux and we ob- 
tain a neutron star temperature of kT°° < 85 eV. By ex- 
trapolating the NSATMOS model fit to an energy range of 
0.01-100 keV, we estimate a thermal bolometric luminosity 

ofi^boi £ lxl ° 33 

erg s . This is comparable to the results 
for the 2003 observation (see Table [3]). 

The 0.5-10 keV unabsorbed flux inferred from the 2009 
spectral data is higher than observed in 2003. This suggests 
that the source intensity may be variable on a timescale of 
years, although the two measurements are consistent within 
the 90% confidence errors (see Table [3]). To further investi- 
gate whether there are intensity or spectral variations be- 
tween 2003 and 2009, we fitted the two data sets simulta- 
neously to an absorbed powerlaw model. When all spectral 
parameters are forced to be the same (i.e., assuming there 
are no spectral differences between the two data sets), we 
obtain N H = (1.3 ± 0.3) x 10 22 cm -2 and r = 1.6 ± 0.3 
for xt = 1-3 (38 dof). Enabling the powerlaw normalisa- 
tion to vary between the two epochs (i.e., allowing only 
the intensity to change) provides a better fit that yields 
N H = (1.2 ± 0.3) x 10 22 cm -2 and r = 1.5 ± 0.3 for 
xl = 1.0 (36 dof). This model fit is shown in Fig. H An 
f-tests suggests that there is a 1% probability of achieving 
this level of improvement by chance. The inferred 0.5-10 
keV unabsorbed fluxes are Fx = (2.0 ± 0.2) x 10" 13 and 
(2.7±0.2) x 10~ 13 erg cm' 2 s' 1 for the 2003 and 2009 data, 
respectively. This suggests that the source intensity changed 
by ~ 30% between the two epochs. 

We next allowed the spectral shape to vary between 
2003 and 2009. The 2000 outburst lightcurye of E XO 1745- 
248 displayed X-ray dips (|Markwardt et alj|2000h . Such be- 
havior has been seen in some LMXBs and is thought to 
result from obscuration of the central X -ray source, e.g., b y 
the outer edge of the accretion disc fe.g.. lBoirin et al.l l2005). 
As a result, the absorption column density along the line of 
sight may be variable. We therefore fitted the quiescent data 
of EXO 1745-248 with the powerlaw spectral index fixed 
between the two epochs, while the hydrogen column density 
was allowed to vary. The results from this fit suggest that 
the X-ray spectrum was possibly more absorbed in 2009, 
although the fit-values obtained for the two epochs are con- 
sistent within the errors (see Table[3]). Keeping the hydrogen 
column density fixed and instead allowing the powerlaw in- 
dex to vary yields a harder spectral index for the 2009 data, 
although within the errors the value is consistent with the 
2003 result (Table [3]). 

The spectral analysis is inconclusive as to whether or 
not the spectral shape changed between the 2003 and 2009 
data sets. We therefore carried out a colour-intensity study. 
For the present purpose we choose two different energy 
bands of 0.5-2.5 and 2.5-10 keV. We determine the ratio 
of Chandra/ ACIS-S counts in the two bands (2.5-10/0.5- 
2.5 keV) and compare this to the intensity in the full en- 
ergy band (0.5-10 keV). The results are displayed in Fig. [5] 
The model-independent colour-intensity study suggests that 
the X-ray emission was harder and brighter in 2009 when 
compared to the 2003 data. For completeness we included 
the 2011 observations (discussed below), but due to the 
large statistical uncertainties we cannot draw any conclu- 
sions about possible changes in the spectral shape compared 
to the other two epochs (Fig. [5]). 
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Figure 3. X-ray count rate lightcurve of EXO 1745-248 obtained 
from the 2009 Chandra observation of Terzan 5, covering the 
energy range of 0.5—10 keV and using a bin time of 1000 s. 
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Figure 4. Chandra/ ACIS spectra of EXO 1745-248 as observed 
in 2003 (bottom, grey) and 2009 (top, black). The solid lines 
represent the best-fit results for an absorbed powerlaw model with 
the hydrogen column density and spectral index fixed between the 
two epochs, while the powerlaw normalisation is variable. 



3.2 The 2011 data: flux estimates and constraints 

Using PWDETECT we obtain a tentative detection with a 
statistical significance of 3.5<r at the position of EXO 1745- 
248 for both 2011 observations. In the 2011-February ob- 
servation there are 11 photons detected within a 1.25" cir- 
cular region centred on the source position. Placing three 
regions of similar dimensions at source-free locations near 
the cluster core suggest that ~ 5 photons are expected from 
the backgr ound. Applying the prescription for low number 
statistics of lGehrela (I1986T ). we infer that 12 net photons are 
detected at 90% confidence level. This implies a count rate 
of ~ 4.0 x 10~ 4 counts s" 1 (0.5-10 keV). 
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In 2011 April there are 13 photons detected at the lo- 
cation of EXO 1745-248, whereas 8 are col l ected from our 
selected background regions. Using iGehrelsl (| 19861 ). we find 
that a total of 15 net photons are detected at 90% confidence 
level, implying a count rate of ~ 3.8 x 10~ 4 counts s" 1 . It 
is possible that the excess of photons is due to a statistical 
fluctuation in the background or a different, overlapping X- 
ray source. In this case, these count rates can be regarded 
as upper limits on the intensity of EXO 1745-248. 

Although the low number of photons does not allow de- 
tailed spectral analysis, we can obtain some estimates on 
the X-ray flux that corresponds to these count rates. We 
group the combined (background-corrected) 2011 spectrum 
to a minimum of one photon per bin and fit it in XSpec us- 
ing W-statistics (an adapted versio n of Cash's statistics that 
allows for background subtraction; IWachter et aT1ll979h . 

Assuming that the observed decrease in intensity was 
completely due to a lower powerlaw normalisation, i.e., 
adopting the average spectral shape inferred from the 2003 
and 2009 data of N H = 1-2 x 10 22 cm -2 and F = 1.5 
(Table E}, we obtain a 0.5-10 keV flux of Fx ~ 1.1 x 
10~ 14 erg cm" 2 s" 1 for the 2011 data (cstat=44 for 33 dof). 
This translates into a luminosity of L q ~ 4 x 10 31 erg s" 1 , 
which would indicate a drop in intensity by a factor of 
~ 20 — 25 compared to the 2003 and 2009 observations. 

Given the hints of spectral variability between the 2003 
and 2009 data sets fSec. 13.1] ). it may be more likely that the 
powerlaw index and/or hydrogen column density was differ- 
ent in 2011. For a realistic range in spectral parameters of 



" 2 and T = 1 
•3.4) x 10" 



- 3, we find a 0.5-10 



14 -2 - 

erg cm s 



for 



N H = (0.5-2.0) xlO 22 cm 
keV flux range of Fx ~ (0.5 - 

the 2011 data, or a luminosity of L q ~ (2 — 12) x 10 3 ' erg s 
We thus conclude that the powerlaw flux varied at least by 
a factor of ~ 6 compared to 2003/2009. If the excess of pho- 
tons seen during the 2011 observations is not related to the 
source, however, these estimates represent upper limits and 
hence the flux variation can be as large as a factor > 50. 

To investigate whether the 2011 data can put further 
constrains on the thermal quiescent emission, we added 
an NSATMOS component to the powerlaw fit for Nh = 
1.2 x 10 22 cm" 2 and V = 1.5. As for the 2003 and 2009 data, 
we assume Mns = 1.4 Mq, Rns = 10 km, D — 5.5 kpc 
and a thermal model normalisation of unity. This way, we 
obtain an upper limit on the neutron star temperature of 
kT°° < 42 eV (cstat=44 for 32 dof). The corresponding 
constraint on the thermal bolometric luminosity, obtained 
by extrapolating the model fit to the 0.01-100 keV energy 
range, is Lq bol ^ 7 x 10 31 erg s" 1 . The estimates for the 
2011 observations are consistent with the constraints ob- 
tained from analysis of the 2003 and 2009 spectral data 
(Sec. CD]). 



4 DISCUSSION 

We use four Chandra observations of the globular cluster 
Terzan 5 to study the quiescent spectral and variability 
properties of the transient neutron star LMXB EXO 1745- 
248, which was active in 2000 and 2011. The observations 
were carried out in 2003, 2009 and 2011. No intervening ac- 
cretion outbursts have been observed from the cluster during 
this epoch. 



LO 

d 



in 



2011 



2009 



2003 



5x10 J 0.01 
ntensity 0.5-10 keV (c s" 1 ) 

Figure 5. Observed Chandra / ACIS-S intensity versus colour of 
EXO 1745-248 during three different quiescent epochs. 



In 2003, EXO 1745-248 is detected at a 0.5-10 keV lu- 
minosity of L q ~ 8 x 10 32 erg s" 1 , while its intensity was 
~ 30% higher 6 years later in 2009 (L, ~ 1 x 10 33 erg s" 1 ). 
Apart from flux variations between the two epochs, the 
source count rate varies up to a factor of ~ 3 during the 
individual observations, i.e., on a time scale of hours. On 
both occasions, the X-ray spectrum is best described by a 
simple absorbed powerlaw model. The source spectrum is 
thus dominated by hard, non-thermal emission and lacks 
the thermal emission component that is seen for many (non- 
pulsating) neutron star LMXBs at such quiescent luminosi- 
ties. This implies that all observed variability can be as- 
cribed to the hard spectral component. Our spectral fits 
suggest that a higher powerlaw normalisation alone can ac- 
count for the increased flux observed in 2009 compared to 
2003. However, a colour-intensity study indicates that the 
X-ray spectrum likely changed between the two epochs, be- 
ing spectrally harder (thus possibly more absorbed) in 2009 
than in 2003. 

While EXO 1745-248 is one of the brightest sources in 
Terzan 5 during the 2003 and 2009 observations, we find 
only a small (3.5cr) excess of photons at the source location 
during the two observations carried out in 2011. For these 
tentative detections we estimate a 0.5-10 keV luminosity of 
L,~4x 10 31 erg s" 1 . If the excess of photons seen during 
the 2011 observations is not related to the source, this should 
be considered an upper limit. We conclude that compared 
to 2003 and 2009, the powerlaw flux had faded at least by a 
factor of ~ 6 in 2011, but possibly as much as a factor > 50. 

By adding a neutron star atmosphere model NSATMOS 
to the powerlaw spectral fit of the 2011 data, we constrain 
the neutron star temperature to kT°° < 42 eV and the ther- 
mal bolometric luminosity to L^ bol < 7xl0 31 erg s" 1 . These 
results are consistent with the constraints inferred from the 
absence of thermal e mission during the 200 3 and 2009 data 
sets (Table [3] see also IWiinands etal] |2005h . 
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Table 3. Results from spectral fitting. 



Year 


N H 
(10 22 cm~ 2 ) 


T fcT e <g x 2 /cstat (dof) F x 

(eV) (10~ 13 erg cm~ 2 s" 1 ) 


L q 
(10 32 erg 


q. bol 

; S -1) 


Thermal 
fraction 


Individual fits: PHABS(POWERLAW) 


2003 
2009 
2011 


1.4 ±0.5 
1.1 ±0.3 
1.2 fix 


1.9 ±0.5 0.2(13) 2 - 2 to!3 
1.3 ±0.3 1.3(22) 2.7 ±0.2 
1.5 fix 44 (33) 0.11 ± 0.04 


8 -°ii'.i 
9.8 ±0.7 
0.4 ±0.2 






Individual fits: PHABS(POWERLAW+NSATMOS) 


2003 
2009 
2011 


1.7±0.7 
1.5 ±0.6 
1.2 fix 


1.9 ±0.5 < 89 0.2 (12) 2.8^^ 
1.3 ±0.3 <85 1.3(21) 3.6 ±1.1 
1.5 fix < 42 44 (32) 0.11 ± 0.05 


lot* 

13 ±4 
0.4 ±0.2 


^ 14 
$ 12 
< 0.7 


< 55% 
<42% 

< 34% 






Simultaneous fit: PHABS(POWERLAW), all parameters 


tied 






2003 
2009 


1.3 ±0.3 


1.6 ±0.3 1.3(38) 2.3 ±0.2 


8.3 ±0.7 






Simultaneous fit: PHABS(POWERLAW), N H and T tied 


2003 
2009 


1.2 ±0.3 


1.5 ±0.3 1.0(37) 2.0 ±0.2 


7.2 ±0.7 
9.8 ±0.7 






Simultaneous fit: PHABS(POWERLAW), T tied 


2003 
2009 


1.1 ±0.3 
1.3 ±0.3 


1.5 ±0.3 1.0(36) 


6.9 ±0.7 
10 ± 1 






Simultaneous fit: PHABS(POWERLAW), N H tied 


2003 
2009 


1.2 ±0.3 


1.7 ±0.3 2.0 ±0.3 
1.4 ±0.3 u.y W 2.8 ±0.2 


7.2 ±1.1 
10 ± 1 







Note. All quoted errors refer to 90% confidence levels. The 2011 data was fitted using W-statistics (an adapted version of 
Cash's statistics that allows for background subtraction), while we employed x 2 -statistics for the 2003 and 2009 data. Fx 
gives the total unabsorbed model flux in the 0.5—10 keV band and L q represents the corresponding luminosity assuming 
D = 5.5 kpc. For the model fits that include a thermal component (NSATMOS) we give the thermal bolometric luminosity 
(L* bol ) and the fractional contribution of the thermal component to the total unabsorbed 0.5—10 keV flux (last column). 



4.1 Constraints on the thermal emission 

Whereas neutron stars cool primarily via neutrino emissions 
from their dense cores and via photon radiation from the 
stellar surface, episodes of mass-accretion cause heating of 
the neutron star. Accretion of matter induces a chain of nu- 
clear reactions that deposit heat deep inside the neutron star 
crust, which is thermally conducted over the entire stellar 
body. In ~ 10 4 yr, a neutron star reaches a thermal steady 
state in which heating due to the accretion of matter is bal- 
anced by cooling via neutrino emission from the stellar core 
and photon radia tion from the surface (|Brown et al.|[l998l ; 
IColpi et al.ll200lh . 

According to the minimal cooling paradigm, there is 
an unescapable rate of neutrino em issions that carries awa y 
energy and cools the neutron star (|Page et al.ll2004l . [2006). 
If a neutron star is particularly massive, however, its large 
central density might open up the threshold for more ef- 
ficient n eutrino emission processes that enhance the cool- 
ing rate (|Yakovlev fc Pethick|[20o3 ). Such neutron stars are 
therefore expected to be relatively cold. In principle, mea- 
surements of the thermal surface radiation can constrain 
the rate of neutrino emissions (i.e., standard/slow versus 
enhanced) from th e neutron star core |Yakovlev et alj|2003l : 
iHeinke et aill2009l ). 



For slow neutrino cooling, deep crustal heating should 
provide an incandescent surface emission that is set by 
the time-averaged mass-accretion rate of the neutron star 
as: L lh hn , = {M ion g )Q„uc/m u (e.g.. IfBrown et all Il998l; 
IColpi et all 1200 ll ). Here, Q nuc ~ 2 MeV is the nuclear 
energy deposited in th e crust per accreted nucleon (e.g., 
lHaensel fc Zdunikll2008l ) . m u — 1.66 x 10 24 g is the atomic 
mass unit and (Mi ong ) is the long-term mass-accretion rate 
of the system averaged over ~ 10 4 yr. If the thermal lu- 
minosity is found to be much fainter, this is indicative of 
a relatively massive neutron star that undergoes enhanced 
neutrino cooling in its interior. We can estimate the thermal 
quiescent luminosity expected from deep crustal heating for 
EXO 1745-248, and compare this with our observational 
constraints. 

IHeinke et all (|2003h estimate a bolometric accretion lu- 
minosity of L acc ~ 3 x 10 37 erg s" 1 by analysing Chan- 
dra and RXTE data obtained during the 2000 outburst 
of EXO 1745-248. For a bolometric accretion luminosity 
given by L acc = (GMns / Rns)(M o1j ) , we can estimate that 
the mass-accretion rate during the 2000 outburst was thus 
{Mob) ~2x 10 17 g s" 1 or ~ 3 x 10" 9 M yr" 1 (assum- 
ing a canonical neutron star with Mns = 1-4 Mq and 
Rns = 10 km). The time-averaged mass-accretion rate of 
the neutron star ((Mi ong )) may be approximated by multi- 
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plying the average outburst accretion rate ((M &)) with the 
duty cycle of the binary, i.e., the ratio of the outburst dura- 
tion and the recurrence time. The recurrence time of EXO 
1745-248 is not well constrained, but the historic activity of 
Terzan 5 provides some rough estimates. 

In the past 30 years 8 distinct outbursts have been ob- 
served from the cluster, of which two could be ascribed 
to EXO 1745-248 with certainty (2000 and 2011; see Ta- 
ble [T}. If this is representative for the long-term accre- 
tion history of the source, the recurrence time of the sys- 
tem would be ~ 11 yr. For a typical outburst duration 
of ~ 2 months (0.17 yr) the implied duty cycle would be 
~ 1%. This yields a time-averaged mass-accretion rate of 
(Mens) ~2x 10 15 g s" 1 (~ 3 x 10~ n M Q yr" 1 ). For these 
estimates, deep crustal heating should generate a quiescent 
bolometric luminosity of Z/ q h bol ~ 4x 10 33 erg s _1 . This value 
would increase if more than two of the Terzan 5 outbursts 
were due to activity of EXO 1745-248. 

The above crude estimate suggests that the quiescent 
thermal luminosity expected to arise from deep crustal heat- 
ing is within a factor of a few of the upper limits in- 
ferred from the 2003 and 2009 Chandra data sets (£ q h b i J$ 
1 x 10 33 erg s" 1 ; Sec. 133) , but a factor > 60 higher' than 
the constraints we obtain from analysis of 2011 quiescent 
data (L q h bo i < 7 x 10 31 erg s" 1 ; Sec. 152)1 . If the low lumi- 
nosity observed in 2011 is intrinsic to the source and not 
caused by temporary obscuration of the X-ray emission (see 
Sec. I4.2|l . the neutron star is thus colder than expected for 
a duty cycle of 1%. This might imply that the source on 
average has a longer recurrence time and that the behavior 
seen in the past 30 years is not representative for its long- 
term activity. When slow neutrino emission processes are 
operating in the core, a thermal bolometric quiescent lumi- 
nosity of iq h boi $ 7 x 10 31 erg s _1 would require a time- 
average mass-accretion rate of (Mong) J$ 4 x 10 13 g s _1 
(< 6 x 10" 13 M yr -1 ). If the observed 2000 outburst is 
typical for the long-term accretion history of EXO 1745- 
248, its duty cycle would need to be < 0.02%, meaning that 
the source should reside > 900 yr in quiescence in between 
subsequent outbursts. 

Alternatively, the neutron star in EXO 1745-248 might 
un dergo enhanced c ore co oling, as was previously proposed 
bv I Wijnands et all (|2005l ). The same conclusion was drawn 
for the other transient neutron star LMXB in Terzan 5, 
IG R J17480-2446, based on its q uiescent thermal luminos- 
ity (|Degenaar fc WiinandslfeoilbT ). This would suggest that 
at least one of the two neutron star LMXBs undergoes en- 
hanced core cooling or that both have a long recurrence 
time. The latter scenario would imply that other historic 
outbursts of Terzan 5 were likely caused by one or more 
additional transient sources. 



4.2 Quiescent variability 

Investigation of the four observations discussed in this work 
reveals that the quiescent emission of EXO 1745-248 varies 
both on short and long timescales. Summarizing: during the 
2003 and 2009 observations the source intensity changed by 
a factor of a few on a timescale of hours. The X-ray flux var- 
ied by ~ 30% between 2003 and 2009, i.e., over a timescale of 
6 years, whereas the powerlaw flux had decreased by a factor 



of > 6 (and depending on the spectral parameters possibly 
as much as a factor > 50) in the 19 months separating the 
2009 and first 2011 observation. There is no evidence for 
variability between the two 2011 observations, which were 
separated by ~ 2 months. It is of note that the source en- 
tered a new accretion outburst f ive months later in 2011 
October (|Altamirano et a l. 2011a). Several (weakly magne- 
tised) transient neutron star LMXBs have been found to 
exhibit changes in their quiescent luminosity and this vari- 
ability has been interpreted in different ways. 

For a group of five sources, KS 1731-260 



JWiinands et all |2002|; ICackett et al] l2010j) MXB 1659- 

l2008h . XTE 



29 (Wijnands et al 



20041: 



Cackett et al 



J170 1-462 llFridriksson et ali |201Q|. |2011|). EXO 0748 



676 (|Degenaar et all 120091 . l201ll : iDfaz Trigo et all l201ll ) 



and IGR J17480-2446 (the other LMXB in Terzan 5; 



iDegenaar et al.ll201ll ) a gradual decrease in thermal X-ray 
emission, extending to many years, is observed after the 
end of an accretion outburst. This variability has been 
ascribed to cooling of the neutron star crust that was 
heate d due to the accretion of matter (e.g., Rutlcdge et al. 
I2002T ). For one of these sources, XTE J1701-462, episodes 
of increased thermal emission were found superimposed on 
the decaying trend. These involved a correlated increase 
in the contribution of the powerlaw spectral component 
and have therefore been attributed to t emporal low-level 
accretion o nto the neutron star surface (|Fridriksson et ahl 
1201(1 1 201 1| ). Neutron star crust cooling cannot be invoked 
as an explanation for the quiescent variability of EXO 
1745-248, since its quiescent emission spectrum is fully 
dominated by hard non-thermal emission. 

Both Aql X-l and Cen X-4 have been exten- 
sively studied in quiescence, using observations spanning 
decades, which has revealed intensity va riations occur- 
ring on time scales of minutes to years |Rutledge et al.l 
19991 . 120021: ICampana fc Stellal 120031 ; ICampana et alj|2004 



Cackett et all l20ld . l201ll ). For Cen X-4, a recent study 



showed that both spectral components need to be variable 
to ex plain the observed intensity variations |Cackett et al.l 
2010). In particular, it was found that the ratio of the power- 
law to thermal flux remained approximately constant, while 
the total flux changed by a factor of > 4. This suggests a 
direct connection between the two spectral components and 
led to the interpretation of residual accretion acting in quies- 
cence <|Cackett et alj20ld ). Studies of Aql X-l were not con- 
clusive about which of the spectral components was underly- 
ing the observed variability. Both accretion onto the neutron 
star surface with a variable mass-inflow rate, as well as the 
variable interaction between the pulsar relativistic wind and 
matter outflowing from the companio n star were suggested 
to explain the observed variability jRutledge et al.l Il999l . 
120011 120021 : ICampana fc Stellall2003l : ICackett et alJl201lD . 

The accreting millisecond X-ray pulsars NGC 6440 X-l 
and SAX J1808. 4-3658 both displayed quiescent flux vari- 
ations by a factor < 2 that could b e entirely attributed 
to the powerlaw spectral component ([Cackett et alj 120051 : 
ICampana et~ai1l2008l ; llieinke et alj|2009h . The behaviour of 
EXO 1745-248 is thus perhaps most reminiscent of these 
two neutron stars, although its variation in powerlaw flux 
appears to be stronger (up to a factor > 6). It is of note 
that some transient LMXBs harbouring a black hole pri- 
mary exhibit X-ray intensity variations in their quiescent 
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state that are of similar magnitude as observ ed for EXO 
1745-248 l|Kong et al.ll2002l ; Irlvnes et alj|2004 ). These sys- 
tems display pure powerlaw spectra in quiescence, without 
any thermal emission component. The mechanism underly- 
ing the variations in powerlaw flux seen in black hole LMXBs 
might also produce the observed variations in the neutron 
star systems. This would necessary imply, however, that the 
powerlaw emission cannot be related to any distinctive prop- 
erties of a neutron star, such as a solid surface and a mag- 
netic field. 

Accretion is known to cause considerable variability on 
long timescales, driving the outburst and quiescence cycles 
of transient LMXBs, as well as on short time scales during 
outburst episodes. Continued low-level accretion therefore 
seems an attractive explanation for the observed variability 
of EXO 1745-248. The source might have been accreting at 
0.5-10 keV luminosities of L q ~ (8-10) x 10 32 erg s" 1 when 
it was observed in 2003 and 2009, while its true quiescent 
luminosity (i.e., the emission level when all accretion activ- 
ity has ceased) is lower and possibly represented by the 2011 
observations (L q ~4x 10 31 erg s _1 , 0.5-10 keV). However, 
it may not be straightforward to explain the large and rel- 
atively rapid drop in X-ray intensity that occurred between 
2009 and 2011 within this framework. This would require a 
strong reduction of the accretion flow in the 19 months sepa- 
rating the two observations, while it had restored again only 
~ 5 months later when the s ource exhibited a new acc retion 
outburst (2011 late-October; [Altamirano et aljfeoilal ). 

Possibly, the non-thermal X-ray emission was not in- 
trinsically reduced but rather obscured during the 2011 
observations. EXO 1745-248 displayed prominent dipping 
activity during the fir st ~ 50 days of its 2000 outburst 
( Mar kwardt et al]|2000l ). During the dips the X-ray intensity 
was reduced by a factor of a few, indicating that the cen- 
tral X-ray source was temporarily obscured, although the 
dips disappeared and did not return in the second half of 
the outburst (Altamirano et al. in preparation). The same 
process may cause obscuration of the X-ray emission in qui- 
escence. For example, the occurrence of a new outburst of 
EXO 1745-248 in late-2011 indicates that there was likely 
a (cold) accretion disk present during the 2011 quiescence 
observations (which were separated by ~ 75 days; see Ta- 
ble^ that may have partly blocked the central X-ray source 
and hence caused the faded X-ray emission compared to the 
2003 and 2009 data. 

Our colour-intensity study and spectral analysis sug- 
gests that the hydrogen column density may have varied 
between the 2003 and 2009 observations, which could lend 
support to the presence of an obscuring medium in quies- 
cence. This interpretation would require the powerlaw emis- 
sion to be produced close to the central X-ray source. An 
origin at the shock front between the pulsar wind and mat- 
ter flowing out from the companion star would then be less 
likely than accretion of matter onto the neutron star's mag- 
netic field or surface. We note that depending on the exact 
spectral parameters, the decrease in quiescent powerlaw flux 
might be as large as a factor of > 50 (Sec. 13.2] ). which would 
be much stronger than the intensity variations seen during 
the dipping activity in the 2000 outburst of EXO 1745-248. 

EXO 1745-248 adds to the (growing) list of neutron star 
LMXBs that show considerable X-ray variability during qui- 
escence. New Chandra observations of Terzan 5 have been 



planned for the coming year. These allow to further study 
the quiescent properties of EXO 1745-248 and can poten- 
tially shed more light on the nature of its strong quiescent 
X-ray variability. 
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